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A First Taste of Dynamical Fermions with an 0(a) Improved Action. 
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We present the first results obtained by the UKQCD Collaboration using a non-perturbatively 0{a) improved 
Wilson quark action with two degenerate dynamical flavours. 
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The Symanzik improvement program is an at- 
tractive method which allows reduction of dis- 
cretization errors order by order in a in physical 
quantities. Recently, the improvement countert- 
erm csw(5o) ^^r the improved action 



S = S'wiison + acsw 



(1) 



has been computed non-perturbatively, thus 
yielding full 0{a) improvement, both in the 
quenched approximation and with Nf — 2 
dynamical fermions . 

The UKQCD Collaboration intends to exploit 
the action (|l|) with two degenerate sea quarks for 
its Dynamical Fermions Project. The calculations 
have been carried out on a Cray T3E recently in- 
stalled in Edinburgh. The version on which these 
results were obtained consisted of 96 processor el- 
ements (PE), each of 900 MFlops peak speed, and 
our code sustained a speed of 25-30 GFlops. 

The algorithm that has been implemented 
is a Generalized Hybrid Monte Carlo (GHMC) 
with even-odd preconditioning. A variant of the 
Sexton- Weingarten integration scheme has been 
employed. We refer the reader to Q for all details 
of the implementation, optimization and testing 
of the algorithm. 

The simulation parameters are summarized in 
tab. 0. The value of csw used was kindly co- 
municated to us by the Alpha Collaboration [|| , 
albeit at a very preliminary stage. Hence the dis- 
crepancy between the number we chose with re- 
spect to the final result presented in ||^. Since 
this is the first attempt to do hadronic physics 
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*Talk presented at Lattice '97, XV International Sympo- 
sium on Lattice Field Theory, Edinburgh (UK), July 1997. 



Table 1 

Simulation parameters. 

with such an improved action in the dynamical 
case, even with slightly incomplete improvement 
we hope to gain much insight in terms of both 
the performance of the code and of the choice of 
the parameters for future runs. Thus, the results 
presented here should be regarded primarily as a 
testing ground to explore the feasibility of the new 
action and algorithm. The choice of a large cou- 
pling was also motivated on such grounds. Given 
the exploratory character of the simulations, we 
have limited ourselved to using local sources and 
have not exploited any form of smearing. Never- 
theless, some light spectrum physics in the mass 
region around the strange quark mass can still be 
extracted. 

For each of the four values of Ksoa we have simu- 
lated valence quarks with the same set of four k's. 
One must keep in mind that in order to take a sen- 
sible chiral limit, we need to put k = K^ca = Hvai, 
but the full (Ksoa, i^vai) plauc can be exploited for 
"strange" physics At /3 = 5.2 the n/p mass 
ratios for the relevant k's are shown in tab. |^. 
These numbers show how relatively far we are 
from the chiral region. 

1. Setting the scale 

The critical point was determined from a lin- 
ear extrapolation in 1/k to vanishing AI^. The 
linear ansatz is well supported by the numerical 
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Figure 1. (a) Lattice spacing and (b) Static quark 
potential as a function of the sea quark mass. 



data, cf. tab. ||. The resuh Kcrit = 0.14033(3) was 
confirmed by the scahng behaviour of the GHMC 
algorithm as a function of Ksea @ ■ 

The lattice spacing was determined in two in- 
dependent ways, cf. fig. from the scale ro 1^ 
and from Mx* 0. We stress that neither method 
requires any extrapolation to /tcrit and are com- 
pletely consistent at all values of Kgca- We view 
the standard approach of fixing the scale from the 
chiral extrapolation of Mp as unreliable, given the 
strong non-linear behaviour in 1/k shown by the 
data, cf. tab. |[ Moreover, in full QCD the p 
can decay and it is therefore dubious to extrap- 
olate to a region where it becomes unstable. In 
fig. we show the static potential from which ro 
was calculated, together with the value of twice 
the static heavy-light meson mass, calculated for 
K = 0.1390, which is the upper limit for the po- 
tential in the case of string breaking. Plotting 
the potential in units of rp, we recover a universal 
curve, as shown in Our data does not exclude 
string breaking although the distances r/ro are 
too small and we need to explore larger distances 
to be able to draw any definite conclusion. 

The strong dependence of a on Ksea indicates 
the non-negligible effect of the improvement coun- 
terterm of the coupling, necessary for the full 
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Table 2 

Some light hadron results for different k's. 



0{a) improvement of the action 
9li'^ + bg{gl)mqa). 



9l 
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To one- loop in perturbation theory bg{gQ) = 
0.012 A^/(?Q 1^, which shows why this countert- 
erm is not needed in the quenched approxima- 
tion {Nf = 0). A non-perturbative method to 
evaluate bg has been proposed in and its ap- 
plication is currently under investigation by us. 
Preliminary results at higher values of /3 indicate 
a much less pronounced mass dependence of the 
lattice spacing, as expected. 

2. Light Spectrum 

We now give some preliminary results of the 
calculation of the light hadron spectrum. The 
standard way to address this point is to show the 
effective mass. In fig. |^ we report it for the pion 
and for the nucleon, for the different k's. The 
lines convey some of the information from the ex- 
ponential fit: the central line is the value of the 
mass, the outer lines denote the error spread, and 
the extension of the lines the fit interval. The er- 
ror has been obtained with the jacknife method. 
The information we can gather from them is that 
the plateau is quite evident for the pion, even 
with a relatively small volume, and the fit is quite 
reliable (a similar situation occurs for the p) , tak- 
ing into consideration the fact that this result is 
obtained without any type of smearing. On the 
other hand, and not surprisingly, the situation is 
worse for the nucleon. In both cases we expect a 
better isolation of the gromid state from smear- 
ing. A full statistics study, including the effects 



of smearing, is underway ||1C | 
suits are summarized in tab 



Some of our re- 
in which we also 



present the variable J = M^'dMp/ dM^. 
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fit time interval= 7-11 (Pion), 7-10 (Nucleoli) 
o a:=0.1370 □ /c=0.1380 * /c=0.1390 ^=0.1395 

Figure 2. Effective mass of tlie pion and nucleon 
for different k's. 

3. Improvement of the Axial Current 



Tlie improved on-shell local axial current 
^ii^, ca) = Af,{x) + acAdf,P{x) 



(3) 



can be determined by enforcing PCAC, as done 
in the Schrodinger Functional formalism An 
equivalent approach, proposed in can be ap- 
plied using standard hadronic correlators. The 
current quark mass (averaged over space) 

i {dtAi{t,cA)pm 

mpcAc(t,c^)=2 (P(t)P(o)) (4) 

is required to be independent of t (up to 0(a^)) if 
the improvement coefficients cgw, ca and bg are 
correctly fixed. We know that this is not com- 
pletely the case for the present simulation, but we 
are here interested mainly in showing the viabil- 
ity of the method. Plotting Atopcac(^*7 c^) = 
mpcAc{T/2,CA) - topcac(**,ca), we fix ca by 
looking at the zero intercept. In fig. ^ we report 
a typical example, which shows that, for a given 
value of t* , there is a distinct cross-over from neg- 
ative to positive values at a fixed value of ca, 
which can be precisely determined by a linear in- 
terpolation. We have checked that other values 
of t* do not yield contradictory information, as 
AmpcAC remains constant with ca within statis- 
tical errors. This analysis suggests a value of ca 
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Figure 3. Atopcac as a function of ca- 

much larger than the perturbative one, which is 
not surprising given the large coupling. Unfortu- 
nately, a comparison with the quenched case is 
not possible since from Q we have ca only down 
to (3 = 6.0 and the (/?, Ksoa) parameters we have 
used are best compared to a quenched value of 
/3 — 5.7, for which the lattice spacing is almost the 
same. A more detailed study, with higher statis- 
tics, is currently underway. Future runs with the 
correct csw and bg are in progress. 

I would like to thank G. Martinelli, G.C. Rossi, 
C.T. Sachrajda, S. Sharpe, M. Testa and all the 
members of the UKQCD Collaboration for fruit- 
ful discussions. We acknowlegde the support of 
PPARC through grant GR/L22744 for the time 
allocation on the Cray T3E. 

REFERENCES 

1. B. Sheikholeslami and R. Wholert, 
Nucl.Phys.B259 (1985) 572. 

2. M. Liischer et ai, Nucl.Phys.B491(1997)323. 

3. K. Jansen, R. Sommer, these proceedings. 

4. Z. Sroczynski et al, these proceedings. 

5. S. Guesken, these proceedings. 

6. R. Sommer, Nucl.Phys.B411 (1994) 839. 

7. C.R. AUton et al, Nucl.Phys.B489(1997)427. 

8. S. Sint and R. Sommer, Nucl.Phys.B465 
(1996) 71. 

9. G. MartineUi et al., |hep-lat/970501^ . 

10. UKQCD Collaboration, in preparation. 



